Buckling of stiff polymer rings in weak spherical confinement 
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Confinement is a versatile and well-established tool to study the properties of polymers either 
to understand biological processes or to develop new nano-biomaterials. We investigate the confor- 
mations of a semiflexible polymer ring in weak spherical confinement imposed by an impenetrable 
shell. We develop an analytic argument for the dominating polymer trajectory depending on poly- 
mer flexibility considering elastic and entropic contributions. Monte Carlo simulations are performed 
to assess polymer ring conformations in probability densities and by the shape measures asphericity 
and nature of asphericity. Comparison of the analytic argument with the mean asphericity and the 
mean nature of asphericity confirm our reasoning to explain polymer ring conformations in the stiff 
regime, where elastic response prevails. 
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INTRODUCTION 



ties are determined by the elastic energy 



It is the interplay of elastic energy and entropy that 
governs the equilibrium form and the dynamics of semi- 
flexible biopolymers. Their competition determines the 
shape and consequently the function of a biopolymer as a 
building block in the cytoskeleton [l|(2j or as an accessible 
storage medium for genetic information [3]. Experimen- 
tal quantification of the elastic and entropic properties 
of biopolymers often employ confinement, may it be by 
clamping one end of the polymer [H [5] or confining the 
whole polymer into a channel [S"-^ or micro-chamber [9] . 
In natural conditions the confinement imposed by cell 
walls and membranes, cell nucleus or viral capsids is ap- 
proximately spherical. This inspired to use the rather 
weak confinement of artificial giant vesicles as a versa- 
tile and well-controllable model system for the investiga- 
tion of polymer and polymer bundle characteristics [10- 
[T2] . Especially but not only in these biomimetic sys- 
tems, that investigate both biological processes and new 
nano-biomaterials, polymer rings become of larger and 
larger importance, stirring theoretical studies of semi- 
flexible polymer rings [13H2Qj . DNA on the one hand 
naturally occurs in ring form [lU [23] while actin and 
actin bundles self-assemble into rings under various con- 
ditions HH [m |22l [111 [25] . Polymer rings are an ideal 
object to investigate entropic and elastic effects as their 
topology induces Euler buckling even in weak confine- 
ment, where the confining cavity is just equal or a little 
larger than the average size of the polymer, see Fig. [ij 
Therefore, spherical confinement serves indeed as an ex- 
cellent tool to investigate the mechanical properties of 
semiflexible polymer rings and how they are affected due 
to biological processes under well-defined conditions. 

Within the wormlike chain model semiflexible poly- 
mers are characterized by their bending elasticity that 
opposes the excitation of undulations from thermal fluc- 
tuations [26 . Representing a polymer of bending mod- 
ulus Kj as a differential space curve y{s) of length L 
parametrized by an arc length s, its statistical proper- 
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where t(s) = dr{s)/ds denotes the tangent vector. The 
competition of elastic against entropic contributions is 
reflected in the material specific persistence length, Ip = 
hz/ki) T, which is just the ratio of elastic bending modulus 
and thermal energy. Comparing this length scale to the 
total length of the polymer gives a measure of polymer 
flexibility L/lp. Polymer flexibility easily varies; there- 
fore, our present study takes it as a variable parameter 
within the region of stiff polymers. It was shown that 
polymer rings due to their topology effectively behave 
about five times stiffer than linear polymers, i.e., their 
stiff regime extends up to ^ 5 [18 rendering poly- 
mer rings advantageous to study elastic responses. 

Semiflexible polymers have previously been the sub- 
jects of investigations under conditions of strong confine- 
ment where the confining cavity is much smaller than the 




FIG. 1: (Color online). Dominant shape of a stiff polymer 
ring without (A) and with (B) spherical confinement. (A) 
Without confinement the first order bending mode excited 
by thermal fluctuations induces a planar ellipse that exceeds 
along its major axis the radius of the corresponding rigid ring. 
(B) Enclosed by spherical confinement the otherwise planar 
ellipse is compressed and Euler buckles into a banana-like 
shape. 
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equilibrium size of the polymer. These conditions arise 
in viral capsids and bacterial envelopes and provoked 
both analytical j27MT] and simulation studies j32H36] . 
While most studies focus on linear polymers viral DNA 
may indeed be circular as taken into account for the in- 
vestigation of knotting probabilities of polymer rings in 
strong confinement [ST . Motivated by nanotechnological 
advances to study polymers in biomimetic systems semi- 
flexible polymers have furthermore been theoretically in- 
vestigated in channels [38^-^ on spherical surfaces [43]- 
f46j and on two-dimensional planes [47l [48] . Concerning 
equilibrium properties it is usually the most likely poly- 
mer conformation that is relevant for biological processes 
and nanotechnological applications. 

Polymer conflguration and form are well- accessible by 
shape parameters based on the radius of gyration tensor 
Q, given by 

Qij = ^ J dsri{s)rj{s) - ^ J dsri{s) J ds'rj{s'). 

(2) 

The eigenvalues and the direction of the eigenvectors 
Ai, i = 1, 2, 3, of the radius of gyration tensor determine 
the spatial extent of a polymer in space. The degree of 
asymmetry, denoted asphericity A, is proportional to the 
normalized variance of the eigenvalues of Q [49 , 
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where A = Ei=i Ai/3 denotes the mean extent. While a 
spherical symmetric object with A^ = A is characterized 
by the minimal value of the asphericity A = 0, a spherical 
asymmetric rod-like object is represented by its maximal 
value A = 1. To measure the degree of prolateness or 
oblateness of an object, the nature of asphericity H is 
defined by [50]: 



4(Ai-A)(A2-A)(A3-A) 



(4) 



The sign of the nature of asphericity is determined by the 
product of the deviations of the eigenvalues from their 
mean and is negative for oblate objects and positive for 
prolate ones. Ranging from S = — ltoZl = l the minimal 
value of the nature of asphericity is attributed to a fully 
oblate object such as a disk, while the maximal one is 
assigned to a fully prolate one such as a rigid rod. 

We use these shape measures to investigate the form of 
stiff polymer rings in weak spherical confinement imposed 
by an impenetrable shell. Employing both Monte Carlo 
simulations and analytical calculations we discern elas- 
tic and entropic contributions and faithfully describe the 
dominant polymer conformation depending on polymer 
fiexibility. In Sec. |ll| we develop an analytic argument 
for the trajectory of the dominant polymer conformation 
considering both entropic and elastic effects. In Sec. Ill 



we assess polymer configurations in spherical confinement 
over ranges of fiexibilities by simulation generated prob- 
ability densities. Finally, we compare asphericity and 
nature of asphericity calculated from our analytic argu- 
ment to their mean values obtained from simulations in 
Sec. 



IV 



In the desired stiff regime our analytic argu- 
ment explains the observed polymer configurations for 
any weak spherical confinement. We conclude in Sec. [Vj 



II. BUCKLING OF AN ELASTIC ELLIPSE 

To understand the form of polymer rings in spherical 
confinement it is insightful to have a description of the 
mean polymer conformation. As the distribution of stiff 
polymer configurations is indeed sharply centered around 
the mean, we develop an analytic argument for the space 
curve of this dominant polymer configuration depending 
on the strength of the confinement and polymer fiexi- 
bility. Based on this dominant space curve (DSC) the 
governing polymer form can be understood and assessed 
by calculating its shape parameters. The successful map- 
ping between DSC and simulation results then also ascer- 
tains our fruitful insights into the whole polymer configu- 
ration. The DSC of a fiuctuating stiff polymer ring arises 
from the interplay of elastic and entropic forces. We an- 
alyze their infiuence subsequently. To derive the DSC of 
a stiff semifiexible polymer ring in weak confinement it 
is instructive to consider first the DSC of an unconfined 
polymer ring. 

A completely rigid polymer ring of contour radius Rc 
is circular. Subjected to thermal fiuctuations, it assumes 
the shape of a planar ellipse [51 , the conformation in- 
duced by the first bending mode. Increasing fiexibility 
enhances the eccentricity of the ellipse within the stiff 
regime. While the major axis of the ellipse grows, the 
minor axis decreases with the square root of the fiexibil- 
ity ^/L/lp p!3HT5] . As spherical symmetry is broken, this 
change in shape yields an increase of entropy and, hence, 
minimizes the free energy. Thus the DSC of the planar 
stiff polymer ring can be parameterized by 



x{s) =Rc I l-^Jf ) 



y{s) =Rc I I + tWt^ ] cos (-|- 




(5) 



z{s) =0, 

where s/Rc G [0, 27r] here represents the polar angle of 
the trajectory and 7 denotes a dimensionless parameter 
that measures the influence of flexibility. The DSC de- 
scribes a polymer ring that is deformed from an oblate 
circle to a more and more eccentric ellipse as the flexibil- 
ity increases. During this growth of eccentricity the to- 
tal length of the space curve is not conserved, hence, the 
model does not predict the overall size of a polymer. This 
caveat does, however, not prevent successful predictions 
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FIG. 2: (Color online). Buckled rod of length S with maxi- 
mal height ho. The distance between its hinged ends 2R de- 
termines the rod trajectory parameterized by i^(s), s G [0,S] 
and the maximal height ho. 



of the shape parameters. As length-invariant measures 
the asphericity and the nature of asphericity are only af- 
fected by the aspect ratio of the axes. In summary, the 
elliptical form of a free polymer is an entropic effect that 
can, however, be translated into an elastic response in 
confinement. 

Confining a stiff polymer ring inside an impenetrable 
sphere induces a change in its sh ape. If the major axis 
of length S = 2Rc{l + ^yfTTjV^) exceeds the diameter 
of the sphere 2i?, no planar ellipse can develop inside a 
sphere. Instead the major axis and therefore the whole 
polymer ring is compressed by the rigid confining walls 
of the sphere into a curved banana-like ellipse as shown 
in Fig. [l] This elastic response of the stiff polymer ring 
to the confinement results in a z-component of the DSC. 
This additional component generates a banana-like poly- 
mer ring, that can be computed by drawing the analogy 
to the buckling of an elastic rod. 

The conformation of a rod of length S pushing against 
rigid walls a distance 2R apart is equivalent to the shape 
of a rod of length S being compressed by a constant force 
/ at its hinged ends. As illustrated in Fig.|2| the configu- 
ration of a buckled rod of length S can be parameterized 
by the angle "di^s) between the tangent vector t(5) and the 
direction r parallel to the compressing force, where the 
arc length s runs from to the length S of the rod. Note, 
that the rod can rotate and move freely perpendicular 
to the axis of the force. Reflecting the mirror symme- 
try of the conformation, the absolute value of the angles 
at both ends of the rod is equal: 'i^(O) = 'i^o = —^{S). 
Euler-Lagrange theory predicts the optimal shape of a 
compressed rod as the state of minimal elastic energy. 
The compressive force / adds a potential term to the 
bending energy of an elastic rod [52^, 
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where we already replaced the bending modulus by its 
relation to the polymer specific persistence length. As 
shown in Ref. the minimization of the above elas- 
tic energy Eq. ([6) under the given constraints results 



in a Euler-Lagrange equation. The minimizing two- 
dimensional space curve describing the optimal filament 
shape is then given by its component r{s) along the direc- 
tion of the force and the component h{s) perpendicular 
to it, see Fig. [2j 
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where a denotes sin('?^o/2) and K, E, and cn are the ellip- 
tic integral of the first and second kind and the Jacobi el- 
liptic function, respectively. When the spatial constraint 
r(0) — r{S) = 2R is respected and the elliptic integrals 
are expanded for small opening angles i^o^ the maximal 
height ho = \h{S/2) — h{0)\ depends on the distance be- 
tween the confining walls 2R and the length of the elastic 
rod 5* only. 
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see Fig. |2j Based on this result the z-component of the 
DSC due to elastic forces can be predicted. Respecting 
the differential continuity of a buckled ellipse the height 
modulation function is taken to be a squared sine result- 
ing in the following z-component for the DSC in confine- 
ment 



z{s) = ho{S,R) sin^ 



(9) 



In addition to the elastic response due to compression 
also entropic forces contribute to the DSC of a spherically 
confined polymer ring. For simplicity, we first assume the 
confining sphere to be of the same radius R = Rc diS the 
contour radius of the confined polymer ring. In this case 
any finite temperature causes the major axis of the ensu- 
ing ellipse to exceed with its apices the spherical confine- 
ment and, hence, forces the polymer to buckle. The elas- 
tic bending energy would be smallest if the ellipse's apices 
both rest on an equatorial plane. Namely, such a config- 
uration maximizes the distance between the apices and, 
hence, minimizes the curvature of the state. Disregard- 
ing rotational symmetry, there is only a single equatorial 
plane. However, entropy increases if the apices may rest 
on any plane instead of just a single equatorial plane. 
This increase in entropy clearly goes at the expense of 
stronger bending. Therefore, the magnitude in deviation 
from the equatorial plane should be related to polymer 
flexibility. As a good estimate we take the DSC of a poly- 
mer ring to nestle half its total height ho{S^ R) below an 
equatorial plane and the other half above. Employing 
Pythagoras la w the total len gth of the major axis is then 
confined to ^/ — {ho/2y. Hence, the ^/-component of 
the DSC follows as 



y{s) = Vi?2-(/io/2)2cos 



(10) 
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Surely, this nestling below the equatorial plane has also 
an effect on the buckling height itself, but it is only of 
second order and, therefore, neglected in the following. 

The above equations for the polymer ring's DSC al- 
ready allow a successful prediction of the mean shape 
of a polymer ring when the confining radius equals the 
contour radius R = Rc. Next our analysis is extended 
to larger radii to enable a full description for any kind 
of confining radius between R = Rc and R = oo. At 
small flexibilities the DSC of a polymer ring is not sup- 
posed to be affected by spherical confinement. As the 
major axis of the ensuing ellipse does not yet stretch be- 
yond the confining walls, a planar ellipse forms as de- 
scribed by Eq. ([5|. Intuitively one would guess that 
confinement effects become noticeable once the length of 
the DSC's major ellipse equals the radius of the sphere 
Rc{l -\- 'y^/L/lp) = R. However, the broadness of the 
distribution of states makes confinement affect the DSC 
even before the major vertices of the DSC's planar el- 
lipse encounter the sphere's shell. So far our considera- 
tions only described the DSC as the mean space curve 
irrespective of the broadness of the distribution of states. 
However, the fraction of configurations with longer ma- 
jor axis forces the DSC to buckle at lower flexibilities 
than expected. This results in an effectively reduced 
radius of the sphere, which we account for by choosing 
R = R — (1 — a){R — Rc). For the lower limiting case 
a = the effective confinement R = Rc instantaneously 
affects the DSC irrespective of the true radius i?, while 
for the upper limiting case a = 1 only the encounter of 
the DSC major axis with the real confinement R = R 
causes an elastic response. Hence, a denotes the per- 
centage of how much below the real confinement radius 
R statistically confinement affects the DSC. a is like 7 
a numerical parameter to be determined from simulation 
data. 

Together these two entropic effects and the elastic 
buckling determine the DSC of a polymer ring of con- 
tour radius Rc in spherical confinement of radius R> Rc. 
For small flexibilities a planar ellipse develops described 
by Eq. ([5|, that is unaffected by the confinement. This 
regime extends up to Rc{l^jy/L/l^) < R-{l-a){R- 
Rc). For larger flexibilities this inequality topples over 
Rc{l + l^LjTp) > R - {l-a){R- Rc) and the DSC is 
described by 

x{s)= i?^(i_^^)sin(5|^), 

^(^) = ^V^-^(i + 7vf)("(i-f)+7v/f) 

X cos (jf;) , 

Based on this analytic argument for the DSC of a stiff 
polymer ring the corresponding shape parameters can be 



calculated and compared to results from Monte Carlo 
simulations. Qualitative accordance with our assump- 
tions for the elastic and entropic forces is gained from pro- 
jections of polymer configurations into two-dimensional 
planes. 



III. 2D PROJECTIONS OF POLYMER 
CONFIGURATIONS 

A discretization of the space curve of the polymer ring 
enables Monte Carlo simulations, which open insights 
into the governing conformations of polymer rings at dif- 
ferent flexibilities. 

To simulate a polymer ring of circumference L the 
Metropolis Monte Carlo method has been employed. The 
polymer ring is modeled as a discrete polygon, that con- 
sists of N segments of fixed length / pointing in the di- 
rection t. The elastic energy of a single conformation 
depends on the direction between successive segments: 
Ti = Nki)T{lp/L) 'ti+i), where the closure of 

the ring is implemented by periodic boundary conditions 
ti = tiv+i. The polymer ring moves through phase space 
by performing crankshaft moves, restricted by the spher- 
ical confinement: Only configurations, which are located 
entirely inside the rigid walls of the sphere, are consid- 
ered for averaging. To collect uncorrelated data, only 
every lO^th of those configurations is considered. We 
sample 10^ configurations for each averaged data point, 
such that the statistical error lies within the ranges of 
the symbols depicted in our graphs. 

To illustrate the form of the polymer rings of radius 
Rc in spheres with R = Rc different flexibilities, the 
probability density of polymer configurations are shown 
in Fig. [3j The position vectors of all samples of polymer 
configurations are mapped on two-dimensional planes 
spanned by two principal axes of the radius of gyration 
tensor in Eq. ([2|, respectively. Ordering the principal 
axes A^, i = 1,2,3, by the magnitude of their corre- 
sponding eigenvalues the largest axis Ai is taken as ref- 
erence axis, and the planes spanned together with the 
intermediate A2 and the smallest axis A3 are considered, 
respectively, to gain insight into the three-dimensional 
configuration space. Considering a planar ellipse that 
buckles due to confinement as discussed in Sec. [n[ the 
plane spanned by the two largest eigenvalues represents 
the planar ellipse and the plane spanned by the smallest 
and the largest principal axis corresponds to the height 
of the buckling polymer relative to the major axis. 

The probability density of polymer configurations in 
the plane spanned by the largest and the intermediate 
principal axis in Figs. (sJA) - [sj^C) reveals the elliptical 
character of the mean shape of the polymer ring. At small 
flexibilities, Figs.jsjA) andjsjB), the polymer trajectories 
are confined to a narrow rim close to the spherical shell 
that broadens with increasing flexibility. With growing 
undulations along the polymer their intermediate axis 
shortens stronger than the larger one. Hence, polymer 
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FIG. 3: (Color online). Probability density and representative snapshots of polymer rings from Monte Carlo simulation data 
for the flexibilities L/lp = 1,3 and 8. As indicated by the cartoons on top, the polymer configurations in the first column 
(A-C) are projected onto the plane spanned by the intermediate A2 and the largest principal axis Ai. In the second column 
polymer configurations (D-F) are projected onto the plane spanned by the smallest A3 and the largest principal axis Ai. The 
gradient in the density of states from high to low is color coded from bright (yellow) to dark (black), the absolute scale of the 
probability density halves starting from 0.0016 onward as the flexibility increases. The snapshots are chosen such that their 
asphericity matches the mean configuration of the observed ensemble. 



configurations resembling an ellipse with higher eccen- 
tricity become more probable. Beyond the stiff regime 
at large flexibilities, Fig. |3|C), the polymer ring exhibits 
compact configurations and looses the character of a pla- 
nar ellipse. In this semiflexible region, the polymer con- 
figurations take a figure-eight shape as indicated by the 
two yellow semicircles in Fig. |3jC). Due to entropic rea- 
sons the eight consists of two circles with different sizes 
for each single polymer configuration [54j, therefore, the 
density distribution is smoothed out in the overlap re- 
gion of the figure-eight. In the flexible regime, the prin- 
cipal axes shrink further with growing flexibility (data 
not shown). However, their ratio remains asymmetric to 
maximize entropy [55 . 

Studying the density distribution in the plane spanned 



by the largest and the smallest principal axis we observe 
bone-shaped probability densities of polymer configura- 
tions, see Figs. [sJD) andjsjE). In these projections the 
density peaks close to the sphere's rim indicate the po- 
sition, where the elliptically shaped polymer configura- 
tions encounter the sphere's shell relative to the equato- 
rial plane. While the largest and intermediate principal 
axis of the elliptically shaped polymer rings in the stiff 
regime map onto the major and minor axis of a buck- 
ling ellipse, the smallest principal axis points toward the 
height of the buckling ellipse. Hence, the width of the 
probability density along the A3 axis in Figs. ^D) and 
[sJe) indicates the maximal buckling height, which in the 
stiff regime is growing with increasing flexibility. An en- 
tirely rigid circular polymer ring would be located in 
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the equatorial plane. With growing flexibility thermal 
fluctuations force the ensuing ellipse to arch out of the 
horizontal equatorial plane, forming a bend. Thereby, 
the major axis of the elliptical polymer ring is clamped 
below or above the equatorial plane. The position of 
the ellipses' apices, which pushes against the conflning 
sphere, is marked by the density peaks in the bone- 
shaped density distribution. The movement of apices' 
positions away from the equatorial plane with increasing 
flexibility is an entropic effect taken into account in our 
analytic argument in Eq. (11). Beyond the stiff regime. 



undulations contract the polymers to a degree, that they 
are no longer forced to undergo Euler buckling but form 
more and more crumpled conflgurations also diminishing 
the polymers' extent along the smallest principal axis. 

The entropic and elastic effects observed in the den- 
sity distributions are in agreement with the analytic ar- 
gument presented in Sec. [llj To substantiate these quali- 
tative observations, the observed shapes of polymer rings 
at different flexibilities in spherical conflnement are quan- 
tifled by the asphericity and the nature of asphericity. 
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IV. SHAPES IN SPHERICAL CONFINEMENT 

The shape of polymer rings is best captured by the 
asphericity and the nature of asphericity as measures of 
the extent of asymmetry and the degree of prolateness 
and oblateness, respectively. Comparing the mean values 
of these shape parameters for free and conflned polymer 
rings displays the dramatic changes in polymer shape due 
to weak conflnement, as shown in Fig. |4j Based on our 
analytic description for the dominant space curve (DSC) 
the shapes of polymer rings are rationalized, and by cal- 
culating exact values for shape parameters of the DSC 
we now show that our analytic argument is in agree- 
ment with the corresponding Monte Carlo data in the 
stiff limit, see Fig. [5] 

Starting from A = 0.25 and E = — 1 for a rigid ring 
the mean asphericity (A) and the mean nature of as- 
phericity (S) of a free polymer ring flrst grow linearly 
with flexibility L/lp in the stiff regime due to the in- 
crease of the eccentricity of the ensuing planar ellipse [51] , 
see Fig. |4] In the semiflexible regime, the free polymer 
evolves into three-dimensional conflgurations and undu- 
lations lead to crumpling that decreases the variance in 
spatial extent. Thereby, the mean asphericity flnally de- 
cays to the exact value for an inflnitely flexible closed 
Gaussian chain of (A) = 0.2464 [56 . In the course of 
this transition the polymer form saturates to a prolate, 
hence, cigar-like shape. In contrast, spherical conflne- 
ment that is small enough to clamp the largest axis of a 
polymer ring provokes the mean asphericity to decline in 
the stiff regime. Only beyond the stiff regime the mean 
asphericity is observed to grow with increasing flexibil- 
ity slowly approaching the value of a free polymer ring. 
Also the linear increase of the mean nature of asphericity 
(E) of a free polymer in the stiff limit is modifled by the 



FIG. 4: (Color online). Monte Carlo simulation data for the 
mean asphericity (A) and the mean nature of asphericity (S) 
versus flexibility L/lp for polymer rings of contour radius Re, 
that are confined by impenetrable spheres of radii R = 1.0 Rc 
to = oo. Relatively weak confinement already induces dra- 
matic changes in the shape of polymer rings. 



conflnement and results in a sigmoidal curve progression 
toward the plateau in the flexible regime. 

The decrease in the mean asphericity for conflned poly- 
mer rings sets in as the ensuing planar ellipse is restricted 
by the conflning shell and buckles into the third dimen- 
sion. As the major axis of the polymer increases with 
flexibility, the buckled polymer conformation gains height 
and, therefore, looses asphericity. This process progresses 
up to flexibilities of L/lp ^ 3. This marks the end of the 
stiff regime deflned by an elastic buckling. The nature of 
asphericity displays that in the stiff regime the cigar-like 
character of the free polymer rings is suppressed by the 
conflnement in favor of more oblate conformations. The 
inflection point of the mean nature of asphericity reflects 
the minimum of the asphericity. Increasing the size of the 
spherical conflnement from R = l.ORc to R = l.SRc re- 
duces the absolute change in asphericity compared to the 
free polymer case. With weaker conflnement the onset of 
the decline of the asphericity is shifted to larger values 
of flexibility, as the ensuing planar ellipse encounters the 
shell only at higher flexibilities. As the distribution of 
polymer extents is broadening with increasing flexibili- 
ties this transition is smoothed out more if the buckling 
sets in at higher flexibilities. Also the character of the 
nature of asphericity changes at the transition, as clearly 
shown by the Monte Carlo data in Fig. [5] If the exten- 
sion of the polymer rings is smaller than the diameter. 
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FIG. 5: (Color online). Comparison of the mean asphericity (A) and mean nature of asphericity (E) versus flexibility L/lp 
calculated from our analytical description in Eq. ( 5|11 ) (light colors) and from Monte Carlo simulation data (dark symbols) 
for polymer rings of contour radius Rc inside spheres of radii R = l.ORc to l.SRc (lower red series). For reference data and 
analytical predictions for a free polymer ring are displayed in each diagram as well (upper blue series). 



the nature of asphericity grows linearly. Its sigmoidal 
character commences at the transition to buckling. 

Beyond the stiff regime, L/lp > 3, undulations start 
contracting the buckled ellipse inducing crumpling to in- 
creasingly compact configurations. Thereby, the polymer 
configurations become less affected by their confinement 
and both shape measures increase toward the value of un- 
confined polymer rings. However, over the range of flex- 
ibilities observed, even the values in the flexible regime 
remain distinct. Although the majority of polymer con- 
formations is coiled up within the sphere, very elongated 
configurations are still discarded and the mean values 
differ from the unconfined case. 

Apart from these qualitative considerations on poly- 
mer shape our analytical predictions for the dominant 
space curve (DSC) in Sec. [ll| can be quantitatively as- 
sessed by comparison to the shape parameters aspheric- 
ity and nature of asphericity. Our predictions for the 
shape parameters depend on two parameters, 7 and a, 
as given by Eqs. ([sj [IT]); in the case of R = Rc only a 
single parameter 7 is needed, as a = by definition. The 
results shown in Fig. [5] are obtained by fitting the param- 
eters 7 [57] and a [58] to both observables for all degrees 
of confinement. Different values for 7 are obtained for 
the asphericity and the nature of asphericity. Such as 
the mean asphericity represents the average shape of all 
polymer configurations such does the fitted parameter 7 
only reflect an average of a whole distribution of param- 
eters. Now both asphericity and nature of asphericity 
have differently shaped, broad and highly skewed distri- 
butions. Therefore, the different results for the fitted 7 
reflect only a range of possible values. However, the fit to 
the stronger peaked asphericity may resemble the average 
growth with flexibility sufficiently well. Altogether, the 
fitted curve for the DSC of polymer rings in spheres with 
radius R = 1.0 Rc and R = 1.3 Rc is in good agreement 
with the simulation results. As our analytical argument 
does not capture the distribution of states, the smooth 



transition from planar to buckled ellipses for R ^ Rc 
shows deviations. There, our argument exaggerates the 
transition in a kink for R = 1.1 Rc and R = 1.2Rc. Con- 
firming the quality of our DSC prediction the analytic 
argument for the nature of asphericity even forecasts its 
sigmoidal character. In the range between R = l.lRc and 
R = l.SRc the transition is again not fully captured due 
to the broad distribution of states, however, the dom- 
inant character of the nature of asphericity is well re- 
flected. Based on our fitted parameters the magnitude 
of all three principal axes of the DSC can be calculated 
for all flexibilities up to L/lp = 3. Having two different 
parameters sets at hand only an estimate of the magni- 
tude is accessible. Recalling that we attribute the fitting 
results for the stronger peaked asphericity a better rep- 
resentation of all possible polymer states, we employ this 
value of 7 to predict for example for R = Rc di maximal 
buckling height at = 3 of about ho / Rc ^ 0.7. If one 
extends the polymer model to account for further micro- 
scopic properties as for instance for torsional stiffness, the 
maximal buckling height is expected to be smaller since 
torsional stiffness increases the elastic energy of out-of- 
plane deformations leaving in-plane bending unaffected. 
Hence, polymers with noticeable torsional stiffness would 
form elliptical shapes due to in-plane modes as observed 
for wormlike chain polymers and thus be compressed by 
confining walls, but the resulting buckling would be to 
less extend. 



V. CONCLUSION 

In summary, the shape and conformation of stiff poly- 
mer rings of contour radius Rc in any weak spherical 
confinement R > Rc imposed by an impenetrable shell 
has been analyzed for varying flexibility L/lp. We find 
that confining a polymer ring induces buckling due to the 
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polymers elastic properties for finite flexibilities L/l^ ^ 3. 
We discern elastic and entropic contributions to the form 
of polymer rings by simulation- derived probability den- 
sities and an analytic argument for the dominating poly- 
mer trajectory. While the elastic response can be sum- 
marized to Euler buckling, the entropic contribution that 
broadens the number of accessible states, induces three 
main effects in the stiff regime. First, entropy promotes 
planar ellipses for any non-zero flexibility, which increase 
in eccentricity with growing flexibility. If eventually the 
major axis is compressed by the confining cavity the poly- 
mer ring buckles as an elastic response. Here entropy 
again takes action as it shifts the plane in which the ring 
is compressed from the energetically favorable equatorial 
plane to smaller radii. At last, due to the broad distri- 
bution of polymer configurations the transition to buck- 
ling is premature and smooth. These four effects are 
sufficient to explain the form of polymer rings in weak 
spherical confinement as shown by comparison of shape 
parameters calculated from our analytic description and 
averaged simulation data. 

Our analytic description, hence, gives a faithful repre- 
sentation of stiff polymer ring conformations and the scal- 



ing of the principal polymer axes, especially the scaling 
of the buckling height, with polymer flexibility. As our 
analytic argument accounts beyond polymer flexibility 
for different radii of confinement, the dominant polymer 
conformation is now available for active control by these 
two experimentally adjustable factors. Employing our 
results in biomimetic experiments, biological processes, 
that are strongly dependent on polymer configuration, 
can be investigated under well defined conditions. Fur- 
thermore, the knowledge of full polymer conformations 
is one of the first steps to build nano-structures based on 
biopolymers. 
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